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O^l , Abstract 

' The non-relativistic scattering of a spin-1/2 particle off a self-dual monopole reduces, 

for large distances, to the dyon problem, studied previously by D'Hoker and Vinet. The S 
matrix (calculated by Zwanziger's algebraic method based on the o(3, 1) © o(3) dynamical 
I symmetry, discovered by D'Hoker and Vinet) is shown to factorize into the product of the 

spinless S'-matrix, Sq, with a spin-dependent factor. The total cross-section is identical to 
the one found in the spinless case. 
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The motion of fermions in a monopole field has been extensively studied [H [21 O HI [5], 
. in particular in connection with the Rubakov effect [1]. For a pure Dirac monopole, the radial 

0^ I Hamiltonian is not self-adjoint, and the solution depends on the boundary condition at r = 0. For 

■ a non-Abelian monopole with symmetry-breaking 'Higgs' potential [6j, the picture is essentially 
Q . the same: the Higgs field approaches its asymptotic value as e~^'^ (where /i is the Higgs mass), 

I and outside the monopole core we are left with a pure gauge monopole, which behaves as an 

O ' Abelian object [6j. 

■ The situation changes dramatically in the Prasad-Sommerfield limit [6] of vanishing Higgs 
potential: the long-range massless Higgs field produces an attractive force (also responsible for 
the existence of static, self-dual, 'Bogomolny' multimonopoles [7j). The test particles can not 

^ ■ reach the origin, curing the non-self-adjointness of the Hamiltonian. The Dirac equation has 

^ , bound state solutions [5]. In particular, there is an interesting E = ground state. For large 

r, the exponential terms can be dropped, the Jackiw-Rebbi equations can be solved exactly; 
electric charge is conserved; the bound state spectrum is of the Coulomb-Kepler type [8]; the 
scattering cross-section factorizes and has a Rutherford form. 

Here we derive similar results in a non-relativistic framework, by enlarging the analogy with 
the Kepler case and, in particular, with its dynamical symmetries [HI'S]. 

We first deduce a generalized Pauli equation. For large distances, the SU(2) symmetry breaks 
down to U(l) yielding a conserved electric charge, and we show that the problem reduces to the 
(at first sight artificial) system considered by D'Hoker and Vinet [3], whose Hamiltonian is 

H = l(p, - eA,f - i + 1^ + ^ - qa,^. (1) 
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Ho = UP^-eA,f-Z + h + —- (3) 



Here the anomalous gyromagnetic ratio 4 arises due to (anti)self-duality, and the potential 
terms appear as remnants of the long-range Higgs field. This system admits a large dynamical 
(super)symmetry, which has been used [3] to derive the bound state spectrum 

^ = ^-^> P=\q\M\ + h,--. (2) 

For p > Igl + 1, the multiplicity of the p-level is 2(p^ — g^), and for the E = ground state, it is 
2\q\. 

Here we focus our attention on the scattering, and calculate the S'-matrix in the spin- 1/2 case 
by generalizing Zwanziger's method [2], based on the dynamical symmetry. We find that the 
spin-1/2 S-matrix factorizes into the product of the spinless S-matrix, So, with a spin-dependent 
factor. 

Zwanziger originally devised his method for the spin-0 counterpart of ([T]), whose Hamiltonian 

is 

1 1 

- + — 1 . 

r 2r2 2g2 

This system has dynamical symmetries analogous to those in the Kepler problem [51 [U] , allowing 
for an algebraic derivation of the bound state spectrum as well as the 5-matrix, 5*0 in Eq. ([7]) 
below, [71|8l|9]. The dynamical symmetry of ([3]) is generated by two conserved vectors, namely 
by the angular momentum, 

V = e'^^x^v^ - q— (4) 
r 

(where = — eA'^), and a "Runge-Lenz" vector, 

A'o = ^e'^'' {v^L^ - Vv'') - ^. (5) 

For fixed value of the energy, E < \ j2q^ , the rescaled Runge-Lenz vector K = q{\. — 2q^E)^^l'^Ko 
extends the manifest o(3) symmetry into a dynamical o(4) algebra. The Pauli-Bargmann [8] 
method yields remarkably the bound state spectrum ([2]) once more, but without the ground 
state -E = 0, and with degeneracy (j? — q^), half of that of the spinning case. 

For energies larger than \/2q^ , the vectors J and K = q{2q^E — l)~^/^Ao close rather to 
o(3,l), with Casimirs 

J.K = M j2_K2 = g2_J__i, (6) 

where k is the wave vector and k = |k|. Then the algebraic method of Zwanziger [2] yields the 
S'-matrix [9], 

So(l,k) = (l(out),k(in)) = 5{E, - E^) ^ {2j + 1) [jj^Hll ^(%,.)(^r'^k)> (7) 



i>l9l 



where E\^ is the energy corresponding to the wave vector k, Ey, = k'^/2 + l/2q^ (so that k = 
(2E — l/q^Y^"^), and ii^ is a certain rotation matrix which brings the z-axis into the direction 
of k. Equation d?]) leads to the modified Rutherford cross-section 



dap 
du) 



4fc4 ) cosec 2^ • (8) 
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In order to study a spin-1/2 particle in a BPS monopole field, let us consider a 5-dimensional 
extended space-time with metric qab = diag((7^jy, 544) = diag(— 1, 1, 1, 1, 1). Let us chose the 
representation 

^ \ Q -il )' ^ V / ' WO 

of the Dirac matrices, so that {7^,7'^} = 2g^^ , A,B = 0,1,2,3,4. Consider a Dirac field ^ 
which transforms according to an unitary representation ^ of the gauge group SU(2), 

and let Qa {a = 1,2, 3, Qa = Qa', [Qa,Qb] = i^abcQc) denote the isospin generators in this 
representation. To a Yang-Mills-Higgs configuration {A^, with the Higgs field $ in the adjoint 
representation we can associate a pure YM configuration Aa in 5 dimensions by setting = 
Afj_, A4 = The 5-dimensional Dirac equation read^ 



^ridf. + a^lQa) + 7'(54 + iAIQa) - mj ^ = 0. (9) 
Introducing 

TT = -iV + QaA", TTi = -idi + Qa^", ^ = 

this equation becomes 

idt(p = {a- ■ TV - iTT4^)x + AQQaf + rrup, (10) 

idtX = (o- • TT + i-Ki) if + A'^QaX - mx- (11) 

In order to take the non-relativistic limit, we introduce, as usual, the slowly varying-in-time 
functions ip and ^, ip = e"*™'*^ ^-iid x = e"™*^. These spinors satisfy 

idtil} = {a ■ -K - 1114)^ + AlQaii, 

idt^ = {(T-7v + i7r4)^ + A^Qa^-2m^. 

In the non-relativistic limit, the second equation is approximately solved by ^ ~ [(T-TT+iTTijip/lm 
(<C V')) so that Eq. ([10]) becomes 

^9t1p = [2k (^^ • ^ - ^^4) (cr • TT + iTT4) + A^Qa] ^. (12) 

This equation can be brought into a more familiar form by developing it as 

(<T • TT — ITT^) (cT • TT + iTT^j = (cT • Tt)^ + i[a ■ TV, 714] + Tt'I, 

and using the identities 

(<T-7r)^ = 7v^ + l[aj,ak][iTj,Trk] = 7v^ + akB^Qa, 

i[cr-7r,7r4] = ak{Dk<^) = ak{Dk<^)''Qa, 

{Bk = \ekijF'^^), which can be verified by a straightforward calculation. Being extensions from 
4 dimensions, none of the fields depend on the extra coordinate, so we require dii) = toci. 
Substitution into Eq. (jl2p yields finally the generalized Pauli equation 



^dti^ = 2k + (Qa^-")' + (yk{Bl + Dk^^)Qa + AlQa] ^. (13) 



^Using this equation means that we choose the Higgs couphng to be the same as the Yang-MiUs coupling, 
more general possibihty would be to require equivariance, —id^ip = C,Tp with some real constant 
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Similarly, in the spin-0 case, one can deduce the Schrodinger equation 



i9t^ = 2k + iQa^")' + ^oQa] i^. (14) 

So far, we worked with an arbitrary Yang-Mills-Higgs configuration in Minkowski space. Let 
us assume henceforth that (A, <I>) is the field of a BPS monopole, i.e., a static, purely-magnetic 
(^0 = 0) solution to the Bogomolny (i.e., anti-self-duality) equations, Bk = D^^, such that 
<I> ^ 1 as r — > 00. Introducing the spin matrix Sk = crk/2, the Pauli equation ([13]) takes the 
form 

i^t^ = 2^ + iQa'^"? + ^^SkB^Qa] i^. (15) 

Notice that the factor 4 in Eq. (fT5]l , the anomalous gyromagnetic ratio of the particle, arises due 
to the Bogomolny equation = D^^. For the spherically symmetric solution of Bogomolny, 
and of Prasad-Sommerfield [6j, in particular, 

= -^(cothr - -), A? = ea^k-{l - -^). (16) 
r r' r smhr' 

For large distances, the exponentially decreasing terms become small and the fields approach 
their asymptotic values, 

^'^ = --(1--), A'i = eaik-. (17) 

The Higgs field is in the adjoint representation, so it can be used to define the electromagnetic 
properties [TU]. Indeed, the electric charge operator, Qem-, can be introduced as, 

Q,^ = ^ = -l=^-Q,. (18) 

All fields are required to be eigenstates of Qem, whose eigenvalues (identified with the electric 
charge) arc integer multiples of ci smallest electric cliarge, [lOj , ^ — '^Qmin-^ ^ — 

0,±1,.... If 

$/|$| is covariantly constant, the electric charge operator commutes with the Hamiltonian and 
the electric charge is conserved. This is what happens asymptotically (i.e., for large r) in the 
BPS case. The gauge potential in Eq. pT|) is in fact a gauge transform of a Dirac monopole 
potential A'^, (rot A )^ = x^/r^. (The above BPS solution has unit magnetic charge g = 1.) 
The electric charge operator is defined hence and is given by 

Qem — Qa' i^^J 

r 

Due to our normalization, qmin = \- Restricted to Qem = gV' states, $'*(5a = q(l — 1/r), B'^Qa = 
—qx^/r^, and thus tt^ = {—idk — qA^)"^. For large distances, the non-Abelian equation ([15]) 
reduces, hence, to 



2m 



-idk - qA^f - — + \ + g2 - 2qak^ 



V'. (20) 



Rescaling as r — > q^r, t — > qH (and for m = 1), this is exactly the D'Hoker-Vinet system ([TJ. 
Similarly, in the spinless case, one gets the Zwanziger system (3). 

Had we studied instead the motion of a spin-1 /2 particle in the field of a self-dual (D<I> = — B) 
monopole, the Pauli term would disappear altogether, and both spin components would satisfy 
independent spinless equations yielding the bound state spectrum ([2]) and the S'-matrix ([7]). 
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In the spin-1/2 case, D'Hoker and Vinet [3] show that the system ([T]) admits three conserved 
vectors, namely 

/ = V + la' (21) 

= Al-e'^k^jyk^la'-qx\^)-{^)a\ (22) 
r r-^ 2q 

= ^a'v^-(^a^vjy + e'i''(^--)a^v''--(^--?a\ (23) 

r q' 2 r q 

The structure involved becomes clearer if we consider 

Ji = r + ^n' (24) 

= {l-2q^H)-^/^{qA' - ^Q'), (25) 

4 = -^n\ (26) 

For a fixed value < E < l/2q'^ of the energy, the operators (I2il)-(l25])-(l26]) span an o(4) © o(3) 
dynamical symmetry algebra, which has been used [3] to derive the bound state spectrum ([2]). 
For scattering off a monopole, the multiplicities are actually doubled, because q and —q are just 
two charge-states of the same particle with non-Abelian structure. It is worth noting that our 
non-relativistic system also exhibits zero-energy bound states found before in the relativistic 
setup [5j. 

Let us now study those scattering states with fixed energies E larger that l/2q^. is now 

= {2q^E - 1)^1/2 -^n'). (27) 
The commutation relations become 

[j(,J^] = ie^^Kjl [J{,K^] = ie^^^K\ [K\K^] = -ie^^Kj^ (28) 

[4,4] = ie,ki4, [4,4] = [JLk']=0, (29) 

showing that Ji and K now span an o(3, 1) algebra, to which J2 adds an independent o(3). 
Using the formulae in Ref. [3j, the Casimir operators are found to be 

Ji-K = M, j2_K2 = g2_ Ji = 3/4, (30) 

with k = {2E — l/g2)i/2. The Casimirs single out the actual representation of o(3, 1) © o(3) 
which operates on the scattering states. The symmetry generators allow for a complete labeling 
of the states, so we have an irreducible representation (spanned, of course, by non-normalizable 
basis vectors). The o(3) of J2 is in its spin-1/2 representation. Representation theory [TT] tells 
us then that among the scattering states there is a convenient "distorted spherical- wave" basis, 
|/c, ji, m, /i), defied by the relations 

H\k,ji,m,fj,) = Ek\k,ji,m,fj,), k > 0; 

{Ji)'^\k,ji,m,fi) = + l)\kji,m,^j.), ji = \q\,\q\ + I, . . . ; 

Jf\k,ji,m,n) = m\k,ji,m,i2), m = -ji, . . . , ji; (31) 

4\k,ji,m,n) = n\kji,m,fi), = 

{k',j[,m',fi'\k,ji,m,fi) = 6{Ek' - E^) 5ji^^j^5„i' ,m^pL' ,1, ■ 



5 



We want to calculate the 5-matrix 5(1, cj|k, s) = (1, (T(out) |k, s(in)). Here |k, s(in)) and |k, s(out)), 
respectively, are solutions of the time-dependent Pauli equation for the Hamiltonian ([T|) , which 
approximate sharply peaked wave packets incoming and outgoing with velocity k at t = ^oo. 
The "helicity" of the incoming (outgoing) particles is fixed to be s = ±1/2, when they are far 
away from the scattering center. This means that, for t = ^oo, the relations 



k, s 
k, s 
k, s 



m 
out 

in 
out 

in 
out 



k, s 



m 
out 



k,s 



m 
out 



(32) 



(cT-k) 



k, s 



m 
out 



2s 



k, s 



m 
out 



with r = r/r, k = k/fe, are verified up to terms which do not contribute to the 5 matrix. This 
implies in turn that our "distorted plane wave-like" scattering states can be taken to satisfy 



H 



k, s 



in 
out 



(Jl-k) 
(K-k) 
(J2 • k) 



k, s 



k, s 



k, s 



m 
out 

in 
out 

in 
out 



Ek 
±q 



k, s 



k, s 



m 
out 

in 
out 



k, s 



k, s 

in 
out 



in 
out 



(33) 



I, a 



m 
out 



k, s 



m 
out 



<5a,.<5(l-k). 



Alternatively, one is (in principle) able to check these formulae on explicit scattering wave 
functions. Anyway, postulating these crucial relations allows us to calculate the S'-matrix. The 
method is to first expand |k, s(in)) (resp. |k, s(out)) ) in the spherical-wave basis ([3T]) and to 
calculate the matrix elements by using these expansions. This is a practical method since the 
spherical waves are easy to handle. Let us choose indeed, for an arbitrary vector k, a rotation 
i?k = exp[ia'^(k)T"], which brings the unit-vector z (directed along the z-axis) into the direction 
of k. Due to spherical symmetry, we can choose scattering states satisfying the phase convention 
relations 



k,s 



in 
out 



exp[ia'^(k)J"] 



kz, s 



m 
out 



(34) 



Here the rotation matrices are replaced by the conserved angular momentum components 
J", which generate the rotations in the quantum mechanical state space. As a consequence of 
the second and fourth Eqs. in ([33|) . |A;z,s(in)) (resp. |A;z,s(out)) can be expanded as 



kz,s i 


in \ 




out / 



ii>kl 



in, out 



k,ji,m = ±q,fj, = s 



(35) 



and all we have to do is to determine the expansion coefficients. (The factors (2ji + 1)^/^ are 
included for convenience.) 
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The third equation in (|33|) tells us how = K • z acts on the l.h.s. of Eq. (|35|). At the 
same time, its action on the r.h.s. of (jSSp is fixed by the representation theory. Putting these 
two pieces of information together, one gets recursion relations for the coefficients b. These are 
actually the same as in the spin-0 case, since the Casimirs of o(3, 1) are identical. We can take, 
therefore, independently of s, the 6's to be those coefficients calculated by Zwanziger [2] for the 
spinless case. 



1 in, out 



The relations J = 
as exp[ia"'{k)J°-] 



Ji + [Ji, J2] = imply that the rotation operator exp[ia"(k) J*^] factorizes 
= exp[ia"(k) Jf ] • exp[ia'^(k) J2 ]. This then gives rise to the equation 



ji 



■m'=-ji fi'=±l/2 



(37) 



with the rotation matrices T> defined by exp[ia"(k) J"] and exp[ia"(k) ] in the representations 
belonging to spin-ji and spin-1/2, respectively. Collecting our results, we obtain the ^-matrix 



5(1, a, k, s) =5{E,- E^)vlil^ {R^'R^) ^ {2j + 1) 

i>l'?l 



(j + i/ky 



^(-,,,)(^r'^k), (38) 



showing that the spin and momentum dependence indeed factorizes. This equation can be 
rewritten in fact as, 

S{l,a\k,s) = V\11^{R-^R^) . 5o(l|k). (39) 

We see that S and 5*0 have the same poles, namely at j — i/k = —n, n = 1, 2, ... , yielding the 
positive-energy bound state spectrum ([2]) with n = n + j = \q\ + 1, . . . . The zero-energy bound 
states do not appear in the S matrix. Equation (I39p implies that, for a particle with incoming 
helicity s and outgoing helicity s', the scattering cross-section is, 

1/2 

where dj ^{9) is the matrix of a rotation by angle 6 around the y-axis in the representation 

spin-1/2. Since s, s' = ±1/2, dJ^2,i/2(^) = ^-V2-i/2(^) = cos ^6*, dY-f/2,i/2(^) = ^i/2-i/2(^) = 
sini6', the total cross section 



—\ =-\^( —\ = ^ (41) 



is identical to that in Eq. ([8]) for spin-0. Remarkably, the results are essentially the same as in 
the relativistic treatment [5]. The factorization of the S'-matrix is actually a general property, 
explained by supersymmetry [12] . 
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